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We outline a run plan and examine the precision with which supersymmetric and 
Higgs parameters may be determined in 1000 fb“^ of running at the nominal 500 
GeV e+e“ linear collider, assuming the Snowmass benchmark mSUGRA point 
SM2 (similar to the SPSl benchmark). 


1 Introduction 

The physics program of the linear e^e“ collider (LC) can be very rich, partic¬ 
ularly in the case that low mass supersymmetry exists. We examine a run plan 
for such a physics-rich case, in which the energies of the LC are selected to first 
measure the sparticle masses approximately using energy distributions at the 
highest available energy, followed by a series of scans at selected superparticle 
pair thresholds. In the supersymmetry part of this study, we pay attention to 
the way in which one can identify specific sparticle states based on the distinc¬ 
tive final states with little background. This paper is an update and extension 
of the study reported in the 2001 Snowmass Workshop proceedings^. 

We assume a physics scenario in which supersymmetry is characterized by 
the minimal supergravity model corresponding to the Snowmass SM2 bench¬ 
mark (mo = 100 GeV, mi /2 = 250 GeV, tan/3 = 10, Aq = 0, and sign^ = +). 
This benchmark point has nearly the same sparticle mass and branching ratios 
as the subsequent SPSla point We take the lightest Higgs boson mass to be 
120 GeV, and to have standard model (SM) properties. 

We have considered the initial 1000 fb“^ of LG running at 500 GeV, ex¬ 
pected to take roughly the first seven years of operation. For runs below 500 
GeV, we assume the luminosity falls as 1 /E. We assume that the LG operates 
with either left- or right-handed electron polarization of 80% and assume no 
positron polarization. The proposed runs and associated polarization states 
are given in Table 1. When a run with left or right polarization is called for, 
we assume that about 90% of the data is taken with the specified polarization 
and about 10% with the other polarization. This run plan differs from that in 
Ref. [1] in that there is no dominantly right polarization run at 410 GeV, and 
there is added running at the Z boson mass. 
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Table 1: Run allocations for the SPSl Minimal Sugra parameters. 


Beams 

Energy 

Pol. 

J Cdt 

[J -^^^jequiv 

Comments 

e^e“ 

500 

L/R 

335 

335 

Sit at top energy for sparticle masses 

e^e“ 

Mz 

L/R 

10 

45 

Calibrate with Z’s 

e^e“ 

270 

L/R 

100 

185 

Scan Xi’ X 2 ’ threshold (L pol.) 

Scan Ti Ti threshold (R pol.) 

e^e“ 

285 

R 

50 

85 

Scan threshold 

e+e“ 

350 

L/R 

40 

60 

Scan it threshold 

Scan CR CL threshold (L & R pol.) 

Scan Xi^ Xi~ threshold (L pol.) 

e+e“ 

410 

L 

60 

75 

Scan T 2 T 2 threshold 

Scan ]l£~ threshold 


580 

L/R 

90 

120 

Sit above xi^ X'^ threshold for X 2 * mass 

e e 

285 

RR 

10 

95 

Scan with e“e“ collisions for cr mass 


2 Initial estimates of sparticle masses from end points 

For the SM2 Susy benchmark, the sparticle masses and primary decay branch¬ 
ing ratios are as given in Table 2. The h^, and masses are 113, 

380, 379 and 388 GeV respectively. The spartners of the lighter quarks and 
the gluino have masses of about 530 and 595 GeV. The two top squark masses 
are 393 and 572 GeV. 

Although many particles are accessible at 500 GeV, the states ti, T 2 , 

Xi^ have dominant decays into r final states. In addition, the sneutrinos have 
dominant wholly invisible decays. On the other hand, the has distinctive 
decays involving Z bosons. Such observations are indicative of the fact that 
the Susy spectra are likely to vary drastically as parameters are changed, and 
thus specific conclusions on the precisions obtainable for sparticle properties 
in one scenario should not be blindly translated to others. 

The traditional ‘end point’ mass determinations seek the sharp edges 
in the energy distribution of a SM particle C in the two body sparticle decay 
A —> BC. Such sharp edges require that a given final state channel under 
observation be fed mainly by just one Susy reaction, and that the SM par¬ 
ticle C be stable. Neither of these conditions hold for many of the channels 
involved in our benchmark Susy point. The second condition is of course not 
actually required, since any known kinematic distribution can be fit with a 
set of templates derived for a set of hypothesized masses, and the best mass 
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Table 2: Masses (in GeV) and dominant branching fractions for the SM2 benchmark. 



M 

Final state 

(BR(%)) 




gr 

143 

Xi"e (100) 





eL 

202 

Xi°e (45) 

Xi^J^e (34) 

X2°e (20) 




143 

XiV (100) 






202 

Xi°M (45) 

Xi^t'M (34) 

X 2 V (20) 



Tl 

135 

XiT (100) 





T2 

206 

Xi°T (49) 

Xr^'r (32) 

X2°^ ( 19 ) 



Ve 

186 

Xl'^^e (85) 

Xi^eT (11) 

(4) 




186 

(85) 

Xi^M^ (11) 

X2°f^M ( 4 ) 



Vt 

185 

X>r (86) 

Xi^rT (10) 

X2°^"t (4) 




96 

stable 





X2° 

175 

T\T (83) 

CRC (8) 

(8) 



X-i 

343 

(59) 

X2°^ (21) 

Xi^Z (12) 

Xi^h (2) 


X4° 

364 

(52) 

Vv (17) 

T2T (3) 

Xi,2-^ (4) 

(6) 

Xr^ 

175 

TiT (97) 

Xim (2) 

xPtu (1.2) 



X2^ 

364 

X2°W (29) 

Xi^^ (24) 

ivi (18) 

Xi^h (15) 

vd (8) 


obtained from the fits. The problem of multiple reactions feeding a particular 
observed channel is more serious, and must be addressed in detail for any given 
Susy benchmark. We comment here on the likely ability to determine sparticle 
masses for the SM2 benchmark. Our estimates of mass precision have been 
guided by previous studies and we scale errors according to the ratio of our 
event samples relative to these studies. 

We restrict attention to channels in which there are only leptons and miss¬ 
ing energy, for which standard model backgrounds should be of little impor¬ 
tance. We have followed the full decay chains of those sparticle pair reactions 
accessible at 500 GeV, weighted by their production cross sections and decay 
branching ratios. The number of events contributed to all 2, 4 and 6 lepton 
final states from all reactions that feed the channel were tabulated for each 
initial electron polarization state. No SM backgrounds are assumed. We have 
looked for the channel/polarization combination in which the production of 
a specific Susy particle are most dominant; in some cases there is no wholly 
dominant channel, and we imagine that a set of coupled-channel studies will be 
needed to get several sparticle masses simultaneously. We discuss the specific 
sparticle template mass determinations in turn. 

The JIr and JIl masses are simple to obtain. The final state Ip with 
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Table 3: Reactions feeding the rr^ final state and the expected number of events before 
acceptance and efficiency losses in 335 fb~^. 


Reaction 

N(eL) 

N(efl:) 

xi^xr 

97,440 

11,229 

Xi°X2° 

29,424 

6,846 

F+Ff 

11,792 

29,547 

r 2 +T 2 (viaxi^r) 

5,716 

2,027 


3,905 

625 


1,395 

428 

r 2 + T 2 (via Xi^^'r) 

1,004 

356 


805 

644 

Xi°X3° 

71 

85 


right e~ polarization has 95% of its contribution from production, and 

the decay JIr decay determines both the Xi^ and ]1r masses. The 

final state final state from left polarized electrons is 95% due to 

production, with one JIl decay to Xi ^nd the other JIr to i^p,Xi^ 

The fj, energy distribution gives MC/Il) in this case. 

The e^e~Ip final state is dominated by the production of ^~l^r 

and eje£. A Snowmass study ^ has shown that through the double difference 
of the distributions for and e~ end points with L- and i?-polarized electrons, 
the masses of cr, cr (and Xi°) can be well determined. 

In benchmark SM2, the final states comingle many production pro¬ 

cesses. Table 3 shows the number of events from various contributing reactions 
for both electron polarizations. The Xi^) X 2 ^ ^-c^cl ti are all large contributors, 
so their masses must be extracted together. There are however some added 
features that help disentangle the reactions. As shown in Fig. 1, for the Xi"*"^!” 
and Tj^rf, the two r’s tend to be back to back, whereas for Xi°X 2 °: the two 
r’s come from the same parent X 2 s-^icl tend to be more collinear. 

For the three dominant reactions feeding the left-polarized TTJf, channel, 
the sparticles involved in specifying the energy distributions are the well- 
measured Xi^jC, and Jl, together with the Xi^)X 2 ° ^cid ti whose masses we 
want to estimate. The energy distributions of the one-prongs from r decays 
are not box-like, but are as shown in Fig. 2. These distributions can be fit 
however for the unknown sparticle masses using templates. For example, with 
the available statistics in 335 fb“^, we find that the statistical error on the ti 
mass is 0.22 GeV when fitting such templates (keeping the Xi^ and X 2 ° masses 
fixed). In a real analysis, all three masses must of course be varied, but more 
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Figure 1: Opening angle of 1-prongs from t decay (a) from reaction Xi^Xi (b) Xi^X 2 ^ 

from left-polarized electrons. 


channels may be employed. 

Despite the fact that no single reaction dominates a particular channel, 
there are several channels that are dominated by the same unknown Xi ^: X 2 ° 
and Ti masses in different proportions. For example the /i in the /irrT channel 
with left polarization arises from JIl ^ 92% of the time (there are 1400 

events before efficiency losses), so gives a reasonable determination of the X 2 ^ 
mass. Other channels involving the same three sparticles are rr (right polar¬ 
ized), fifiTT (left), fiTTT (left) and tttt (left). A careful study is needed to 
see how well the three unknown masses can be disentangled from a coupled 
channel analysis. I guess that precisions of about 1 GeV should be achievable 
from the 500 GeV runs, sufficient to locate the energy at which a scan should 
be performed. 

The higher mass gaugino masses may be partly measurable. The x^ has 
the useful and distinguishable decay modes, Xs^ Xi°(X 2 °)-^ from which the 
X 3 ° mass can be determined, albeit with limited statistics. The x^ mass 
though accessible kinematically, has event rates that are too small for a mass 
measurement. The X 2 ^Xi^ reaction threshold opens at 539 GeV in the SM2 
benchmark; if we can deduce that this is is likely to be the case, we would 
raise the LG energy to about 580 GeV to allow X 2 ^Xi^ production and give a 
rough estimate of the X 2 ^ mass. 

The z7 masses are difficult to obtain in this scenario. The most favorable 
channel for isolating Vg pair production is the errr final state, where VgV* 
contributes about 40% of the 6.5K events (before efficiency losses). The eerr 
and e/i/ir channels have significant VgV* contributions. Perhaps it is possible 
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Figure 2: Energy distributions of 1-prongs from r decay in the t+t channel; (a) all 1- 
prongs, and (b) the more energetic 1-prong. 


to dig the Vf, mass out. The and Vr seem hopeless to measure, as these 
sparticles do not come close to dominating any final state channel. 

3 Sparticle threshold scans 

Previous studies ^ have examined the precision that can be attained on mass 
measurements from the appropriate energy scans near a two particle thresh¬ 
old. For this measurement, the problem of disentangling the reaction feeding 
a particular final state does not occur, providing that there are no overlapping 
thresholds in a particular channel. However, cross-sections near thresholds 
tend to be small, so the scans can use rather significant portions of the lu¬ 
minosity budget, and should be chosen for only those cases where there are 
substantial gains to be made. 

The first studies ^ of scans used ten energy points across the threshold in 
question. Subsequent studies have questioned this strategy, finding that 
particularly in the case of low uxBR and relatively slow threshold turn-ons 
characteristic of sfermion pairs (/3^), fewer points on the scan may be better. 
An analytic study ^ that did not include the effects of backgrounds showed 
that the strategy for minimizing the errors shows little improvement for more 
than about three points on the scan. 

One can improve the slope of the threshold curve for enCR by producing 
them in s—wave from an e“e“ initial state (rather than in p—wave as for e“'"e“) 
even after the inclusion of beamsstrahlung effects and the reduced luminosity 
in e“e“ operation^. We have adopted this strategy in our run plan. 
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Table 4: Mass precisions for the run plan of Table 1. 


sparticle 

SM 

end point 

SM 

scan 

SM 

combined 

Sr 

0.19 

0.02 

0.02 

bl 

0.27 

0.30 

0.20 

V'R 

0.08 

0.13 

0.07 

ML 

0.70 

0.76 

0.51 

n 

~ 1-2 

0.64 

0.64 

T2 

- 

1.1 

1.1 

Ve 

~ 1 

- 

~ 1 


7? 

- 

7? 


- 

- 

- 


0.07 

- 

0.07 

X-i 

~ 1-2 

0.12 

0.12 

X3° 

8.5 

- 

8.5 

X4° 

- 

- 

- 

Xi^ 

~ 1 - 2 

0.18 

0.18 

X2^ 

4 

- 

4 


The estimated precisions were scaled from those of previous threshold scan 
studies. We have not explicitly included the effects arising from differences in 
background at our benchmark point relative to previous studies. However, we 
conservatively considered only the dominant channel and polarization state. 

4 Susy mass and parameter precisions attainable 

The mass precisions obtained for the run plan of Table 1 are shown in Table 4. 
Given these mass precisions and making the assumption that the supersymme¬ 
try is indeed mSUGRA, the errors that can be obtained for the Susy parameters 
are shown in Table 5. 

5 Higgs and top quark studies 

The Higgs production processes occur at most of the energies in our run plan; 
there as many ZH events over the full run as would be obtained in about 
525 fb“^of running at 350 GeV, or in about 1220 fb“^at 500 GeV. We do not 
expect that running at several energies should materially affect the Higgs pa¬ 
rameter determinations, and we simply scale the precision obtained in previous 
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Table 5: Errors on mSUGRA parameters. 


parameter 

SPSl 

mo 

100 ± 0.08 GeV 

Wi/2 

250 ± 0.20 GeV 

^0 

0 ± 13 GeV 

tan/3 

10 0 0.47 


studies by the statistics of the ZH sample. We have used the somewhat 
more favorable estimates for thge TESLA study®, but note that there is still 
some uncertainty on the branching ratio errors, in particular for cc decays. The 
projected Higgs parameter errors are shown in Table 6. Use of Higgs bosons 
produced through WW fusion may improve the situation somewhat. 

The top quark parameters are taken from the scan near the ti threshold 
at 350 GeV. The statistical errors are expected to be small compared with 
the theoretical errors associated with QCD theory. We may expect that the 
mass and width of the top quark may be determined to 150 MeV and about 70 
MeV respectively in the 40 fb“^allocated to this threshold. (Less luminosity 
would be needed to saturate expected theory errors for the top quark, but the 
it threshold overlaps those of 6^6^ and for which larger luminosity 

accumulations are desired.) 

6 Calorimeter calibrations 

It will be necessary to make good relative calibrations of the calorimeter energy 
scale in a LC detector, so as to capitalize upon the excellent intrinsic resolution 
of the calorimeter. These should be done with a periodicity that is shorter than 
the gain-drift time for the calorimeter. Two methods could be envisioned: the 
first employing special runs at the Z boson mass, and the second using data 
collected during ordinary running at higher energies. 

Our run plan in Table 1 calls for four runs of 2 fb“^each at the Z (every 
other year). Others would ask for more - perhaps one per year. For each such 
run, we expect 3.6 million Z ^ ee or and 76 million Z hadrons. We 
assume a calorimeter with of order 10® EM towers and 5 x 10^ hadronic towers, 
and suppose that we calibrate 2x2 blocks of these as fundamental units. We 
take the calorimeter energy resolutions to be 6E/E = O.lb/y/E 0 0.005 (EM) 
and 6E/E = OA/y/E 0 0.01 (hadronic). Then we obtain about 28 electrons 
in each EM calibration block and about 24K hadrons in each hadronic block. 
In each block, the precision on the Z mass measurement is given by SM = 




Table 6: Relative errors (in %) on Higgs mass, cross-section, total width, branching ratios 
and Yukawa couplings (A) for the run plan of Table 1. 


Parameter 

error 

Parameter 

error 

Mass 

0.03 

Ttot 

7 

a{ZH) 

3 

^ZZH 

1 

a{WW) 

3 

5.WWH 

1 

BR{bb) 

2 

\bH 

2 

BR(cc) 

8 

^ccH 

4 

1 

+ 

5 

^ttH 

2 

BR(gff) 

5 

5~uh 

30 


cm/VN, with anf/M 5E/E ai E = Mz- The resulting mass measurements 
for each block yield a statistical error on the Z mass, and hence a calibration 
accuracy of 0.38% (EM) and 0.03% (hadronic). Since these precisions are 
small compared with the intrinsic energy resolution smearing, they should be 
adequate for measurements of SM and Susy masses. 

It is highly desirable to have calibration methods that can be performed 
during high energy running, and thus can track changes in calibration. Bhabha 
scattering has been studied as a means of continuous EM calibration. Here 
we focus on an alternate method using ‘energy flow’ (not to be confused with 
the technique of substituting charged particle momentum measurements for 
calorimeter energy deposits). We define energy flow to be the relative number 
of events seen above some arbitrary energy threshold, over a ring of cells at 
fixed azimuth. The azimuthal symmetry for physics processes means that 
observed variations are due to differences in the relative gain constants ai 
for the set of cells in the ring (Etrue = aEmeas). A simplified estimate of 
the calibration precision obtainable can be made assuming that the energy 
spectrum in a given azimuthal ring is given by dN/dE = Ae~^^. Then for N 
events above an energy threshold Eth, one can show that 5a = l/{BEth\/N), 
and the optimum (smallest 5a) choice of threshold is Eth = 2/3. 

For 10 fb“^of data at 500 GeV and the inclusive particle spectrum (we 
actually used the flatter 1000 GeV Monte Carlo spectrum and hence con¬ 
servative for our purposes), and with 200 cells at fixed (p for each of 20 bins 
in cos 9, we find that a 1% calibration can be obtained. Though of poorer 
precision than the special Z pole calibration, this is a very useful check for 
drifts in the calibration during high energy running. The method itself can be 
checked with energy flow measurements at the Z pole, where we estimate a 
precision of about 0.8% for each cell. Cross calibration of different azimuthal 
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Table 7: Profile by year of the luminosity accumulation. The luminosity is given in fb ^ as¬ 
suming 500 GeV operation. 


Year 1 2 3 4 5 6 7 

i Cdt 10 40 100 150 200 250 250 


rings would be performed using Z bosons from the high energy runs. 

The absolute energy calibration of the calorimeter at the 0.03% level is 
desired to match the expected statistical precision of the Higgs boson mass 
measurement. This can be obtained after relative calibration of the calorimeter 
blocks in the special Z pole runs discussed above, and should yield calibration 
of both EM and hadronic calorimeters to the requisite accuracy. 

7 Run order 

We should note that before the LC run plan is constructed, we will have 
an understanding of the actual physics scenario from the LHC/Tevatron. In 
particular, the Higgs mass will be known and the presence or absence of su¬ 
persymmetry will be established. The run order given here assuming the SM2 
Susy scenario discussed above is not well tuned, but is offered to guide our 
thinking. Clearly the initial stages will also be constrained by the progress of 
the accelerator complex. We assume the growth of luminosity in the LC to be 
that given in Table 7. 

1. In the first year, if 500 GeV operation is feasible, it would be profitable 
to run there to get a first indication of the masses of the easier Susy 
particles to observe. In the 10 fb”^, we expect ~700 ZH events at 500 
GeV. If 500 GeV running is not available, operation at 350 GeV would 
be called for (~1500 ZH events). If it is only possible to run at low 
energy, then one would devote the initial running to the Z pole. 

2. Go to 500 GeV as soon as possible for ~80 fb“^to obtain about two 

times the ultimate errors on sparticle masses. This should be sufficient 
to establish the eL,eR,Ji,L,]lR, and Xi masses well, and give reasonable 
estimates for ^nd ti masses. 

3. Scans at 285 GeV threshold) and at 350 GeV {tt,eLeR, and 

Xi^Xi~ thresholds). 

4. Gomplete the 500 GeV run for the ultimate precision on end point masses, 
and for disentangling the multi-r final states. 
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5. Scan at 270 GeV {X 1 X 2 ^ 1^1 thresholds). 

6 . Scan at 285 GeV using e~e~ initial state to give the precise cr mass. 

7. Operate at above the nominal top energy at 580 Gev with e+e“to do the 
Xi^X 2 ^ study and obtain the X 2 ^ mass. 

Note that no Giga-Z run is called for our SM2 Susy scenario, since the 
need to provide the extra Z-pole (and more precise W mass) is not so great in 
the scenarios where supersymmetry provides the mechanism for EWSB. 

8 Considerations for a run scenario if it’s not Susy 

If the world we inhabit does not include low mass supersymmetry, we will know 
this at the start of the LG operation. In this world, there are two avenues for 
exploring the character of EWSB: observe new states or deviations in cross- 
sections from SM predictions at high energy; and improve our knowledge of 
the precision electroweak observables at the Z pole and WW threshold. The 
run plan in this case should reflect these two paths, and we outline such a plan 
below. Approximate equivalent luminosity acquisitions (the accumulations 
that would have occured in the same time at 500 GeV) are given in parentheses. 
We have not specified the order in which the runs might be taken, but clearly 
those involving very careful understanding of the beam energies (for the WW 
threshold run), positron polarization (for Giga-Z), or 77 collisions should be 
left for later in the program when the machine is well-understood. 

1. (500 fb“^) Look for new high mass phenomena 

2. (150 fb~^) Run at the highest accessible energy by trading off luminosity 
for energy. Seek new signals at higher energies and new physics energy 
dependences. 

3. (100 fb“^) Obtain the full Giga-Z sample, with polarized e+. About 20 
fb“^are needed at the Z (thus £equiv ~ 100 fb“^). 

4. (150 fb“^) Run with 77 collisions at the maximum energy available to 
give an alternate view to e“'"e“collisions for the trilinear gauge couplings, 
heavy Z', or complementary sensitivity to large extra dimensions. 

5. (80 fb“^) Scan the WW threshold for a precision W boson mass mea¬ 
surement. 

6 . (20 fb“^) Scan the tt threshold for precision top quark mass and width 
information. 
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9 Conclusions 


In the delightful case that Nature provides us with low mass supersymmetry, 
the LC should offer the possibility for extensive and precise measurements to 
point the way for understanding EWSB and the character of physics beyond 
the SM. The details of the run plan will depend critically upon the exact model 
that we encounter. Considerable ingenuity may be required to unscramble the 
effects of many accessible supersymmetric particles, but in the r-rich case 
examined here, it seems possible to do a good job with this. Even for the 
low mass Susy situation studied here, much running at 500 GeV is needed to 
establish a good picture of the various sparticle masses. And it remains clear 
that further studies at energies above 500 GeV will be needed in this, or other 
conceivable physics scenarios. 
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